Glassy Spin Dynamics in Non-Fermi-Liquid UCu5_a,Pda;, x = 1.0 and 1.5 
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Local /-electron spin dynamics in the non-Fermi-liquid heavy-fermion alloys UCus-j^Pdj,, x — 1.0 
and 1.5, have been studied using muon spin-lattice relaxation. The sample-averaged asymme- 
try function G{t) indicates strongly inhomogeneous spin fluctuations, and exhibits the scaling 
G{t,H) = G{t/H'<) expected from glassy dynamics. At 0.05 K 7(a;=1.0) = 0.35±0.1, but 
7(a;=1.5) = 0.7±0.1. This is in contrast to inelastic neutron scattering results, which yield 
7 = 0.33 for both concentrations. There is no sign of static magnetism > 10^^ fJ-s/V ion in either 
material above 0.05 K. Our results strongly suggest that both alloys are quantum spin glasses. 

PACS numbers: 71.27.-fa, 75.30. Mb, 76.60. Cq. 



Magnetic resonance |]-|| and other |J| experiments 
have demonstrated the importance of structural disor- 
der in the breakdown of Landau's Fermi-hquid theory 
in certain /-electron intermetallic compounds and alloys. 
Disorder-driven mechanisms have been considered for the 
non-Fermi-liquid (NFL) properties of some of these sys- 
tems and it is natural to consider the possibihty 
of extremely disordered or "glassy" behavior. On theo- 
retical and experimental grounds it is known that glassy 
dynamics lead to long-time correlations with distinct sig- 
natures as the freezing or "glass" temperature Tg is ap- 
proached from above In a spin glass the spin au- 
tocorrelation function q{t) = {Si{t)-Si{0)) is theoreti- 
cally predicted to exhibit power-law {q{t) = ct~") or 
"stretched-exponential" {q{t) — cexp[~{Xt)^]) behav- 
ior |10|. Power-law correlation has been found in spin- 



glass AgMn using muon spin relaxation (/iSR) ^ . 

This Letter describes evidence from ^SR experiments 
that spin correlations in the NFL alloys UCu5_a;Pda;, 
X = 1.0 and 1.5, are indicative of glassy spin dy- 
namics. The sample-averaged muon relaxation func- 
tion (asymmetry) G{t, H) is strongly sub-exponential, 
indicating a quenched inhomogeneous distribution of re- 
laxation rates, and obeys the time-field scaling rela- 
tion G{t,H) = G{t/H'') for applied magnetic field H 
between ~15 Oe and ^^1 kOe. The field dependence 
corresponds to a measurement of the Fourier transform 
of q{t) over the frequency range ^^H/2ti k, 200 kHz- 
14 MHz, where 7^ = 27r x 13.55 kHz/Oe is the muon 
gyromagnetic ratio. Power-law behavior of q{t) is im- 
plied by the observation 7 < 1 and also by the 
temperature-frequency scaling found in the inelastic neu- 
tron scattering (INS) cross section [|ll|, although the pos- 
sible connection with glassy dynamics has to our knowl- 



edge not been previously noted. The present measure- 
ments extend by three orders of magnitude the frequency 
range over which power-law correlations are observed in 
UCu5_2:Pda;. Zero-field /iSR above 0.05 K shows no sign 
of static magnetism or spin freezing in UCu4Pd; this to- 
gether with the glassy scaling points to some form of 
quantum spin-glass behavior [ p^Jl3| . 

In /iSR experiments spin-polarized positive muons are 
implanted into the sample, and the subsequent time evo- 
lution of the muon polarization is monitored by detecting 
the asymmetric distribution of positrons from the muon 
decay [Q. Muon relaxation in a magnetic field applied 
parallel to the muon spin direction (longitudinal field) is 
due to thermally-excited /-electron spin fiuctuations that 
couple to the muons. A muon at a given site experiences 
a time- varying local field Hioc(i) due to fluctuations of 
neighboring / moments. Following Keren et al. |^, under 
motionally narrowed conditions the local muon asymme- 
try G{t, H) relaxes exponentially: 



G{t, H) = exp [-2AVc(i7)t] 



(1) 



where — 7^(|HiocP) is the time-averaged mean- 
square coupling constant in frequency units, and the local 
correlation time Tc{H) is given by 



rc{H) 



dt q{t) cos{j^Ht) = cUc{H) . 



(2) 



We consider A and the prefactor c but not the functional 
form of Uc{H) to vary from site to site in a disordered 
material. Then the sample-averaged asymmetry G{t, H) 
is given by 



Git,H) 



dAdcp{A,c)exp[-2A'^cuciH)t] , (3) 
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FIG. 1. Field dependence of sample-averaged muon asym- 
metry relaxation function G{t) in UCus.sPdi.s, T = 0.05 K. 
Curves: fits as described in text. 



FIG. 2. Data of Fig. |l| plotted versus the scaling vari- 
able t/H'^-\ 



where p(A, c) is the joint distribution function for A and 
c. It can be seen that the field and time dependence 
enter Eq. (||) in the combination Uc{H)t. This means 
that G{t, H) scales as this combination independently of 
the form of p(A, c). 

For both the power-law and stretched-exponential 
forms of q{t) this time-field scaling results in 



G{t,H) = G{t/H^) 



(4) 



after Fourier transforming q(t). Here 7 = 1 — a < 1 for 
power-law correlation, and ^ — 1 + f3 > 1 for stretched- 
exponential correlation as long as the muon Larmor fre- 
quency LJfj, = 7^iJ is much greater than A Q . If Eq. is 
obeyed a plot of G{t, H) versus t/H'' will be universal for 
the correct choice of 7. The sign of 7 — 1 distinguishes be- 
tween power-law and stretched-exponential correlations. 

Samples of UCu4Pd and UCus.sPdi.s were prepared 
as described previously jl^. Arc- melted ingots were 
crushed into powder under acetone. The powder was 
passed through a 90 micron sieve, and pressed with GE 
7031 varnish into pellets 13 mm dia. x '^1 mm thick- 
ness. /iSR data were taken at the LTF facility of the 
Paul Scherrer Institute, Villigen, Switzerland, for tem- 
peratures between 0.05 and 1.1 K and for magnetic fields 
between zero and 10 kOe applied in the direction of the 
muon spin polarization. 

Figure| shows G(i) in UCua.sPdi.s for T = 0.05 K and 
values of apphed field H between 13 Oe and 2.5 kOe. The 
relaxation slows with increasing field. For low enough 
fields we expect the field dependence to be due to the 
change of ui^ rather than an effect of field on q{t); a break- 
down of scaling would occur for high fields where this 
ceases to be true. The same asymmetry data are plotted 
in Fig. 1^ as a function of the scaling variable t/H'^ . For 
7 = 0.7 ± 0.1 the data scale well over more than three 
orders of magnitude in t/H'^ and for all fields except 



2.5 kOe. Fields fisH > ksT would be expected to affect 
the spin dynamics, and indeed the static susceptibility 
of UCu4Pd is suppressed by fields '^l kOe below ~0.5 
K (VoUmer et ai, Ref. ||l^). The scaling exponent 7 
is less than 1, implying that q{t) is well approximated 
by a power law (or a cutoff power law [^) rather than 
a stretched-exponential or exponential. From our data 
q{t) w ct^'^-^^'^-^ . We note again that no specific form 
for the muon asymmetry function has been assumed. 

A scaling plot is given in Fig. ^ for UCu4Pd, T = 
0.05 K. Here the scaling exponent 7 — 0.35 ± 0.1 is 
significantly smaller than in UCua.sPdi.s. Scaling again 
breaks down for high enough fields; the data for 2 kOe 
clearly fall above the low-field scaling curve. Data taken 
in UCu4Pd at T = 0.5 K (not shown) scale with the same 
exponent. This is in contrast to the scaling behavior 
of spin-glass Ago. 995 Muq. 005 ||], where 7 varies strongly 
with temperature as the glass temperature Tg — 2.95 K 
is approached from above. 

The fluctuation-dissipation theorem |l^ relates Tc{H) 
to the imaginary component x"{lli) of the local (q- inde- 
pendent) /-electron dynamic susceptibility: 



rc{H) 



(5) 



for hu! <^ ksT. INS experiments |Q show scaling of 
the sample-averaged x"(a;, T) as uj~'^F{hLj/kBT), with 
7 = 0.33 and F{x) = tanh(a;/1.2) for both UCu4Pd 
and UCu3.5Pdi.5. Using this form Tc{H) obtained from 
Eq. (||) is independent of T and proportional to the 
latter is in accord with the /xSR scaling. The INS value 
of 7 agrees with /itSR data for UCu4Pd (7 = 0.35) but 
not for UCu3.5Pdi.5 (7 = 0.7), suggesting that in the lat- 
ter sample the behavior of Tc changes in the unexplored 
region between the highest muon frequencies (~10 MHz) 
and the lowest INS frequencies (~10 GHz). Clearly high- 
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FIG. 3. Scaling plot of G{t,H) vs t///"-^^ for UCu4Pd, 
T = 0.05 K. 
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resolution INS studies (neutron spin echo or backscatter- 
ing) would be desirable between these frequencies. 

To go further one must fit the fiSK data to an appropri- 
ate functional form for the asymmetry. We have chosen 
the stretched-exponential 



G(t) =exp[-(Ai) 



(6) 



where K < 1 gives sub-exponential relaxation corre- 
sponding to a distribution of relaxation rates. This func- 
tion is purely empirical. It is used because it charac- 
terizes an a priori unknown relaxation-rate distribution, 
and because the rate A conforms with a general defini- 
tion of a characteristic rate by the time 1/A where G{t) 
decays to 1/e of its initial value. Equation (||) fits the 
data in applied field to within the statistical error. 

For H — the data were fit to the product of Eq. (j^) 
and the zero-field Kubo-Toyabe (K-T) function char- 
acteristic of static relaxation by nuclear dipolar fields at 
muon sites. This form is expected when the muon lo- 
cal field has both static nuclear dipolar and dynamic /- 
moment contributions. A nuclear dipolar field ~2.3 Oe 
was measured in both alloys for T ^ 1 K, where the 
contribution of U-moment fluctuations to the zero-field 
muon relaxation rate vanishes. Nonzero values of H 
were chosen large enough to "decouple" the muon relax- 
ation from the nuclear dipolar field (l8| leaving only the 
dynamic U-moment contribution, so relaxation data for 
these fields were fit to Eq. (||) without the K-T function. 
Curves giving these fits are plotted in Fig. |l|. 

Figure | gives A(T) and K{T) for UCu4Pd at three 
values of H. Below 1.1 K A increases slowly and satu- 
rates to a constant below 0.1-0.2 K. As noted previously 
a temperature dependence of A(T) is not in agreement 
with the temperature-independent relaxation predicted 
from INS scaling. The exponent K is approximately 0.7 
at 0.05 K, indicative of a broad distribution of relaxation 



FIG. 4. Temperature dependence of muon stretched-expo- 
nential relaxation rate A (filled circles) and exponent K (open 
circles) in UCu4Pd. (a) Applied longitudinal field H — 0. 
(b) H = 16 Oe. (c) H = 300 Oe. 



rates and decreases slightly with increasing temper- 
ature. Similar behavior is exhibited by A(T) and K{T) 
in UCu3.5Pdi.5 (data not shown), with rates slower than 
in UCu4Pd by ~30% at low fields and -100% at 100- 
300 Oe due to the larger scaling exponent. 

No anomaly is found in the zero-field data at temper- 
atures ~0. 1-0.2 K, where specific heat and ac suscepti- 
bility (in different samples) suggest spin-glass-like freez- 
ing ]l6[ |. The muon-/-moment coupling in UCus-^Pdj; 
is predominantly dipolar with a coupling field 0.55 ± 
0.05 hOe/fiB- Randomly-frozen moments of the or- 
der of 1 hb/V ion would result in a muon relaxation 
rate — 50 /is~^, two orders of magnitude larger than 
the observed rate. This result places an upper bound of 
'^10"'^ /^s/U ion on any frozen moment in UCu4Pd or 
UCus.sPdi.s. We believe that the discrepancy with the 
results of Refs. [|l^ results from differences in annealing 
conditions; fiSH experiments to explore this question are 
currently underway. It is noteworthy that the saturation 
of A and K occurs in the same temperature range as the 
reported spin-glass-like freezing |p^ . 

Our results can be compared with existing theories of 
disorder-driven NFL behavior. Preliminary analysis has 
indicated that the order of magnitude of the experimen- 
tal rates at these very low temperatures cannot be di- 
rectly accounted for by the simple single-ion Kondo dis- 
order model of NFL behavior |1£ ] 
of the "Griffiths-phase" theory ^|, 
feet of /-moment clustering, there is no dissipation in the 
/-electron spin dynamics, and the local cluster dynamic 
susceptibility is sharply resonant at a distributed charac- 



In an early version 
which treats the ef- 
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teristic tunneling energy E: 

x"(w, E) oc S{uj - E) tanh(£;/2T) . 



Together with the distribution function P{E) oc 
where A < 1 is a nonuniversal scahng exponent, this im- 
mediately yields 



j dEP{E)x"{uJ,E) (XLu 



-l+A 



tanh(w/2r) 



in agreement with the INS data for A « 0.7 and reminis- 
cent of our /iSR results. 

But the observed time-field scaling of G{t, H) demon- 
strates that it is the local ^"{u}) itself, not just the aver- 
age which scales as uj~^ [cf. Eqs. This is 
not a property of the Griffiths-phase model, in which the 
scaling is found only after the average has been taken, 
and only for a sharply resonant (and nonscaling) form of 
x"(a;, E). A recent form of this theory |20j considers dis- 
sipative effects, which broaden x"(u;,i?) but do not give 
it a scaling form. Furthermore, if the width of x"(o-') is 
much greater than uj it is not hard to show that x"(w) 
no longer follows P{u)—E) so that this mechanism for 
scaling of x"(w) is lost. It is also difficult to see how 
the dynamic susceptibility of Ref. would yield the 
observed temperature dependence of A(aj^,T). Thus the 
Griffiths-phase theory does not seem to account for our 
results for a number of reasons. 

Muon relaxation in UGu5_a;Pda; clearly indicates 
glassy dynamics, and the increase of A(T) at low temper- 
atures (Fig. suggests a low-temperature critical point. 
The increase in K(T) with decreasing temperature might 
be due to more efficient averaging over the disorder as- 
sociated with a growing correlation length. But A(T) 
saturates below ~0.2 K and there is no evidence for a 
true phase transition. The relaxation rates are broadly 
distributed {K < 0.7); we are never dealing with critical 
behavior in a homogeneous system. The distribution of 
rates is itself static, at least over the time scale of the 
experiment (~10 /xs), otherwise G{t) would tend to be 
averaged to an exponential with the average rate. 

The faster relaxation (i.e., slower spin fluctuations) for 
X = 1.0 than for a; = 1.5 implies that the relevant critical 
point is near the concentration for which the Neel tem- 
perature Tn is suppressed to zero {x w 1). But in spin- 
glass AgMn 7^1 from below as T ^ Tg whereas 
here j{x=1.5) = 0.7, 7(1.0) = 0.35 < 7(1.5). A mean- 
field model of a disordered Kondo alloy at a quantum 
critical point predicts j — l/2atT = 0, suggestive 
of our results but not in detailed agreement with them. 

By definition NFL behavior is a property of the lowest- 
lying excitations of a metal, to which a low-frequency 
probe such as /iSR is extremely sensitive. Many aspects 
of the spin dynamics in UCu5_a;Pda; are poorly under- 



stood, but the strong disorder, glassy behavior, and ab- 
sence of a phase transition strongly suggest that these 
(7) alloys are quantum spin glasses. 
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